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Abstract 

We discuss how the mass of new physics particles involved in a pair of short decay 
chains leading to two invisible particles, for example slepton pair production, followed by 
the decay into two leptons and two neutralinos, may be measured in central exclusive pro- 
duction (CEP) with forward proton tagging. We show how the existing mass measurement 
strategies in CEP may be improved by making full use of the mass-shell constraints, and 
demonstrate that, with around 30 signal events, the masses of the slepton and neutralino 
can be measured with an accuracy of a few GeV. 

1 Introduction 

Many theories beyond the Standard Model (SM) include a dark matter (DM) candidate as part 
of the mass spectrum. Typically, the stability of the DM candidate is the result of an assumed 
global Z2 parity, under which the new (SM) particles are odd (even). For example, in the 
minimal supersymmetric extension of the SM (MSSM), this is the so-called R-parity T^p, while 
in models with large extra dimensions, this could be a Kaluza-Klein (KK) parity. 

If such a DM candidate is produced at the Large Hadron Collider (LHC), its presence could be 
revealed through the observation of large missing transverse energy. The parity implies that 
the DM candidate will be pair-produced at the LHC. This leads to complications in measuring 
its properties, in particular its mass, since only the vector sum of the transverse momenta of the 
two DM candidates in an event can be measured. 

Strategies to circumvent this problem abound in the literature (ll-[23| Q Many of these assume 
the existence of certain long decay chains. An important example in the MSSM is the pair 

^ See also Ref. [25| for a review. 



production of squarks (g), followed by the decay chain 

^ /^X? , (1) 

where X12 neutralinos (x? is the DM candidate), and / is a slepton. The existence of 

more constraints, e.g. mass-shell conditions and dilepton invariant mass edge measurements, 
than undetermined 4-moment8[§ would allow simultaneous determination of the mass of the new 
particles. These long decay chains imply various challenges for the relevant mass measurement 
techniques, e.g. correct identification of the decay topology and combinatoric problems in as- 
signing the identified final state particles to appropriate stages of the decay chains, and so it is 
advantageous to consider alternative methods to complement existing strategies. 

The central exclusive production (CEP) of a system X can be written in the form 

pp ^ p + X + p, (2) 

where the '+' signs represent the presence of large rapidity gaps. In this type of reaction, 
only the system X is produced in the central detector, with no additional hadronic activity 



in the absence of pile-up, see for example Ref. (26|, |27| (and references therein) for a review. 
The outgoing protons scatter at small angles and remain intact, so that additional kinematic 
information can be obtained by installing proton tagging detectors far from the interaction 



point [28[. In particular, in the presence of these detectors the momenta of the outgoing protons 
can be measured precisely, and therefore the longitudinal momentum and invariant mass of the 
central system X can be accurately determined. 

In this paper, we present a mass measurement method utilising the additional longitudinal 
kinematic information in the exclusive two-photon production of a pair of charged particles 
{X^,X2 ), 

-> p + 77 + p , (3) 

77 ^ xtx^ , 

followed by the decay 

X± ^ x^ + x, (4) 

where x is a visible particle (or more generally, fully reconstructed system of particles), and x 
is an invisible particle that is not seen in the detector. For concreteness, we will consider the 
specific case of the exclusive two-photon production of a pair of sleptons {If, 12)^ followed by 
the decay into leptons {li,l2) ^^e lightest neutralinos (xi;X2) (note the different notations 
compared with Eq. ([T])), i.e. 

pp p + 77 + p , 

77 lp2 , 

If!^ ltl2XlX2 , (5) 



^In the present context, there arc 6 (= 8 — 2) unknowns: the four momenta of the two x\ (8), minus the 
measured missing transverse momenta (2). 
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leading to a signature with two tagged protons, an opposite-sign same-flavour lepton pair and 
missing 4-momentum with no additional jet activity. The irreducible SM background comes 
from the production of W pairs, i.e. 

pp p + 77 + p , 

77 ^ W+W,- , 

W^W^ -> . (6) 

Compared to the long decay chain discussed in Eq. ([T]), the two-photon production process 
benefits from a clean, purely leptonic signature, with no combinatoric ambiguity. Furthermore, 
assuming the slepton only decays into the lightest neutralino, the process is essentially a two 
scale problem that depends on the (true) masses of the slepton and the lightest neutralino, 
^^tme^ ^toue-j^ together with discrete parameters such as the particles' charges and spins. 

The irreducible SM background is also under good experimental control. This is because 
in the case of the signal, two-photon production must create two same-flavour sleptons. In 
the absence of large lepton flavour violation, this will lead to two same-flavour leptons in the 
final state. For the SM background there is no such constraint on the flavour of the final-state 
leptons: a measurement of the opposite flavour final states e^/i^ would therefore allow a cross 
check of the background determination. 

However, the cross section for the exclusive two-photon production of slepton pairs is small 
when compared to the inclusive production cross section from hard inelastic scattering, which 



is roughly two orders of magnitude larger [29|, |30|: as a result, mass measurement in this central 



exclusive channel is only realistic at the 14 TeV LHC for mj up to about 200 GeV. However if 
the coloured sparticles, i.e. the squarks and the gluinos, are too heavy to be produced copiously 
at the LHC, the presence and properties of light colour singlet sparticles may be best measured 
via two-photon production. If (as we hope) the coloured sparticles are produced in abundance, 
and long decay chains are observed, the present method might still be used as an independent 
check or to provide external constraints on the long decay chains. 

The system in Eq. ([5]) has four unknowns: only the vector sum of the 4-momenta of Xi and 
X2 is measurable. Given a consistent mass hypothesi^ (m^-, m-j^), the four quadratic mass shell 
constraints allow the system of equations to be solved analytically, up to a two-fold degeneracy. 
As a result, a (m|-, m^) mass region consistent with the measured 4-momenta can be obtained on 
an event-by-event basis. We believe this consistent mass regions contains all information that 
can be derived from the measured 4-momenta. Neglecting finite detector resolution effects, the 
(different) consistent mass region from each event must contain the true mass point (m^''™, m|''"°). 
Combining the consistent regions from a number of events will develop a peaking structure 
around this point, hence in principle allowing determination of the true masses. 

In determining the true masses and in separating the signal from SM background, we find 
the maximum values of (m™*^) and my_ (m™*") within the consistent mass region particularly 
useful. This is because, by construction, m™^^ and m™^^ are bounded from below by m^™" and 
^toue respectively. They are also bounded from above by half the 77 invariant mass {m^^) and 



^The technical definition of a consistent mass hypothesis wiU be given in the next section. 
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Figure 1: Cross sections for pair production of spin-1 (y^V~), spin-1/2 and spin-0 

{S^S~) charged particles in central exclusive two-photon production as a function of the particle 
mass (m'''"'^) at the 14 TeV LHC, calculated using the PhoCEP Monte Carlo event generator. 



half the missing energy {E„ 
likelihood functions from m? 



An additional advantage is the relative ease in constructing 
, which could be particularly useful given a small data 
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sample, as is likely the case in light of the small two-photon production cross section. 

The main aim of this paper is to study how well the slepton and neutralino masses may be 
determined using mf^^^ and m^'" computed in each event. These masses can also be estimated 
from m^^ and i^miss |3l|, 12^ , which have threshold values at 2mj^^° and 2m|'"° respectively. We 
shall compare the two methods and illustrate the improvements that can be made by making 
use of the mass-shell constraints. The potential for measuring properties of the SM Higgs 
boson HQ, long-lived gluinosjpj, stops jsej, MSSM Higgs (alH, NMSSM Higgs jssl and 
other BSM Higgs scenarios in the central exclusive mode with tagged forward protons 

have also been discussed in the literature. 

Note that the decay chain considered can be applied to new physics models with different 
particle spins, for example the pair production of charged spin-1/2 fermions followed by the 
decay into neutral spin-1 bosons, as might be expected from a KK model. In fact, it can be 
shown that, given the same masses and charges, the production cross section for spin-0 bosons 
is the smallest when compared with spin-1/2 fermions and spin-1 bosons, see Fig. [H In this 
paper, we shall concentrate on the first case, which is the most conservative, and use the SUSY 
decay chain in Eq. as our working example. 

The paper is organised as follows. In Section [51 we present the analytic solutions consistent 
with the 4-momenta measurable in Eq. and derive the quantities m?'^^ and m^^^. We then 
present the numerical results in Section[3l and compare between the performance of (m™^^, rn^^^) 
and (m^^, i^jniss) in determining (m^'^"°, m^^"*^), after including the CEP dynamics simulated by 



a new event generator PhoCEP [41[]. This is followed by an estimate of 



using the 
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likelihood method advertised above assuming an integrated luminosity of 300 fb at the 14 TeV 
LHC, which is the main result of our paper. We conclude in Section |H 



2 The consistent solutions 

The 4-momenta of the particles in the decay chain in Eq. ([5]) are subjected to the constraints 

PU = PXr+Ph^ 

pI = 

pi = ^1' c^) 

where [i = 1, 2) and p^^ is the 4-momentum vector sum of the two 7's. For simplicity, the 
sign of the sleptons and leptons are dropped. In the above expressions, only p^^ and pi. are 
directly measurable. The 4-momenta of the SUSY particles, pi and Py.., can be expressed as 
linear combinations of p^^, pi- and P, a space-like 4-momentum vector defined by 



which has the property 



Pf, = ef.^paP'^^^pf.pl , (8) 



P-p,, = P-pi, = P-pi, = 0. (9) 



Eq. (171) implies 



Ill 

Px2 = ■^^'^ + a)Py^ + -{b - l)pi, - -{c + l)pi,, + dP , 
Ph = - a)Pii - - l)pi, + ^{c - l)pi^ - dP , 
PI, = l{l + a)p^^ + ^{b-l)pi,-^{c-l)pi, + dP, (10) 
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with mass-shell conditions 

111 

--(1 - a){b + l)p^^-pi, + -(1 -a){c-l)p^^-pi^ - -{b+l){c-l)pi^ ■ pi^ , 

ml, = i(l + a)X7 + ^'^' 

+ ^(1 + a){b- l)p^^ ■ pi^ - ^(1 + a)(c + l)p^^ ■ pi, - ^{b - l)(c + l)pi, ■ pi, , 

ml = l(l-a)X. + rf^^^ 

-^(1 -a){b- l)p^^ ■ pi, + ^(1 - a)(c - l)p^^ ■ pi, - ^(6 - l)(c - \)pi^ ■ pi, , 

m? = i(l + a)X7 + ^'^' 

+ ^(1 + - 1)^77 ■ Vh - ^(1 + a)(c - \)p^^ ■Vh-\^- l)(c - ■ Vh . (H) 

where the parameters a, 6, c and d are to be determined. For simplicity, we have neglected the 
lepton mass mi throughout, although we note that consistent solution regions can also readily 
be defined for general mass values. 

Here m^ and m-^ need not be the true masses, but input parameters consistent with a positive 
c?^. The structure of the system is such that upon subtracting the different expressions in Eq. (fTTI) . 
the additive d^P^ term and the terms bilinear or quadratic in (a, 6, c) drop out. This leads to 
three equations linear in (a,6, c). Defining the quantities 



Am^^ Ai = ^^^-^, A^^^Ii:^, K,,^-^, (12) 

Vh ■ Ph Ph ■ Ph Ph ■ Pl2 Ph ■ Ph 



where 

Am?- = mj-m^, (13) 
the solutions for the parameters (a, b, c) are given by 

(A2 - Ai) 



^^A^^ - 2A1A2 
1 - A2 + Aa 



A^^ - A2(Ai + A2) 
A^^ - 2A1A2 



Inserting Eq. into Eq. ( ITTj) . all four mass constraints give 

A^ = CaAl + CbAA + Cc + d^Ap, (15) 
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where 

Ph ■ Ph Ph ■ Ph Ph ■ Ph 

and 



m% m? p2 
Ax = ^, ^, Ap = , (16) 



1 (Ai + A2)' - 2A^^ 
~ 4 A^^-2AiA2 ' 

Cb = -(A1 + A2-2), 

Cc = ^(A^-,-2AiA2). (17) 

Note that in Eq. f lT5|) . c^, c;,, Cc and Ap are constructed from directly measureable 4-momenta. 
A mass hypothesis (m^-, m^) is consistent with the measured 4-momenta if it results in a positive 
dp, in which case a two-fold degenerate solution for {p^-^^,p^^,pj_^,pjj corresponding to ±\d\ is 
obtained. 

The boundary of consistent solutions on the {mj, m^) plane is given by Eq. f lTS]) by setting 
d = 0. To see this, note that since P is space-like, Ap < 0. This means that for a given Aa, A^ 
decreases as d"^ increases. An upper boundary on the (Am?_,m|) plane therefore has c?^ = 0. 
The consistent region is bounded from above by this quadratic curve, and below by the x-axis, 
i.e. m| = 0, which can be transformed into a boundary in the {mj, ruj^ plane. Within the 
boundary, d'^ > 0. By construction, the solutions are time-like. 

More observations on the properties of the consistent solutions can be made. First, the sign 
of the energy component of the two solutions (for ±|(i|) must be the same, since it is always 
possible to boost to a frame where the energy component of the space-like vector P is zero, 
in which case the two solutions have the same energies. Second, since the consistent solutions 
are continuous functions of and mj, it follows that the energies of all consistent solutions 
must have the same sign. The energies must then be positive, due to the fact that the true 
masses must be a consistent solution, and hence lies within the boundary. Third, the solution 
is clearly symmetric when exchanging Ai and A2. This means that the boundary obtained for 
the mass-shell conditions for p^^^ and p^^j must be the same, consistent with the fact that both 
constraints lead to Eq. (fT5|) . 

As discussed in Section [T], the consistent mass region may be characterised by the maximum 
values of and mj, i.e. m^^^ and m^^^. They correspond to the stationary values of A^ and 
A;- as functions of Aa in Eq. f lTSj) . The results are 



\2 



(mrf = iPk -Pi.) X U - ^^^^ ) . (18) 



By construction, m™*™ and m™^^ are bounded from below by the true and mj respectively. 

So far we have shown that analytic boundaries containing mass hypotheses consistent with 
the measured 4-momenta can be obtained on an event-by-event basis. In the next section, we 
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Before cuts (fb) 


After cuts (fb) 


SM 
SUSY 


1.00 
0.15 


0.72 
0.13 



Table 1: Cross sections for the SUSY model with (mf"",m|™^) = (150, 100) GeV and the SM 
background at the 14 TeV LHC before and after the lepton cuts —2.5 < rji < 2.5 and pt > 10 
GeV. Cuts on the fractional momentum loss of the outgoing protons 0.0015 < .^j < 0.15, are 
applied in both cases. Branching fraction of the W into one lepton flavour is included. 

present examples showing the densities of consistent solutions given a large number of events, 
and show the extent to which the quantities m™*^^ and m™*^^ can be used to determine the true 
masses, crucially in a small event sample. 



3 Numerical study 



Our numerical study is performed using the new Monte Carlo (MC) event generator PhoCEPj4l| . 
It simulates the two-photon exclusive production of a W^W~ pair and a charged scalar pair, 
S^S~, via the decay chains given in Eqs. ([5]) and (Q, as well as the production of two charged 
fermions, F~^F~, although we do not discuss this process in detail in this paper. The equivalent 



photon approximation (28|, |42|] is used to calculate the photon luminosity dL^^/dW due to 
photon emission from the incoming protons, and this is then combined with the known 77 — X 
(X = W^W~, S^S~, F+F~) subprocess cross sections to give the full CEP cross section 

dapp^pxp _ f da^^^x{W)dL^l^_^ ^^^^ 



dfl J dn dW 

where W is the 77 cms energy. The produced particles are then decayed according to the decay 
chains given in Eqs. ([5]) and ([6]), and in the case of W^W~ production full spin correlations for 
the subsequent decay are included. The MC generates unweighted events, and full kinematic 
cuts on the final state particles and outgoing protons can readily be implemented. The survival 
factor 5*^, representing the probability that the outgoing protons do not undergo soft rescattering 
and therefore spoil the exclusivity of the final state, is estimated to be roughly 90% for two- 



photon exclusive production [27[, and for simplicity is not included in the MC. The differential 
cross sections for the 77 — )■ X subprocesses are summarised in Appendix |Al 

In the following, we present numerical results of the methods discussed in Section [2j For 
concreteness, we focus on a SUSY model 

• SUSY model: (mp^m|"°) = (150,100) GeV 

and compare with results for SM W^W~ pair production, which constitutes an irreducible 
background to the SUSY process. For the SUSY model, the slepton is assumed to always decay 
into a lepton and a neutralino. The branching fraction of the SM W into one lepton flavour is 



taken as 10.8% 43 



In our simulation, a final state proton is assumed to be tagged with 100% efficiency if the 
fractional momentum loss lies within the range 
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300 




I ' ' ' 1 

50 100 150 200 

[GeV] 

Figure 2: Density of the consistent mass regions for a SUSY (solid red) and SM (dashed blue) 
sample. From the outermost to the innermost, each set of contours delineate parameter regions 
with probability density larger than 0.1, 0.3, 0.5, 0.7 and 0.9. The mass parameters for the 
SUSY model are (m^''™, ml™") = (150, 100) GeV. No detector smearing is included. 



• 0.0015 < < 0.15 z = 1,2, 

which corresponds to the expected coverage under the assumption that roman pot detectors are 
installed at a distance of 220 and 420m from interaction point, see for example Ref. 44 . 
The following lepton pseudo-rapidity (rji) and transverse momentum (px) cuts 

• -2.5 <rii< 2.5, 

• Pt > 10 GeV , 

are applied to the opposite-sign same-flavour lepton pairs. We assume 100% detection efficiency 
if a lepton passes the above cuts. The cross sections before and after the lepton cuts at the 14 
TeV LHC are displayed in Table [H 

In Fig. m the density of the consistent (m^-, m^) regions are shown. In both the SUSY and 
SM samples, a peaking structure located in the vicinity of the true mass point develops. Finite 
detector resolution effects will smear the spike structure slightly. However, given the LHC's 
ability to measure accurately the 4-momenta of the forward protons (in the presence of proton 
taggers) and the charged leptons, the smearing effects should be relatively mild. The position 
of the high density regions are significantly different such that the SUSY signal can be easily 
separated from the background. 

Next we show in Fig. [3] scatter plots for m^^^ and mf^^^, each of which contains 10,000 events. 
Recall that, by construction these two quantities are bounded from below by m|™° and m^^^. 
As we can see, in both the SUSY and SM samples (mp*^^, m™^^) cluster close to the true mass 
point {mf^^,m^™^). This indicates that a good estimation of the true mass parameters and 
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Figure 3: Scatter plots showing the (mp^,m™^^) values consistent with each event. Left: SUSY 
(^tme^^toue^ = (150,100) GcV. Right: SM. Each plot contains 10,000 events. No detector 
smearing is included. 



identification of the presence of a SUSY signal from the SM background could in principle be 
achievable. 

To show the enhancement obtained by utilising the additional kinematic information, in 
particular in mass determination, we compare our results with an alternative method which 



measures the end points of the missing energy (-Emiss) and m^^ distributions |3l|- The end point 
values of these distributions are at 2m|'"'' and 2m^~^^^ respectively. These distributions take into 
account the sum of longitudinal momenta and energy of the neutralinos that can be measured 
in CEP with tagged forward protons, but do not include information from mass constraints. 

In Fig. m we show the m™*^^ and m™^^ distributions. The -Emiss and m^^ distributions are 
displayed in Fig. O Compared to the m™^^ and m^^^ distributions, the -Emiss and m^^ distri- 
butions are generally broader. This qualitative difference is due to the fact that both _Emiss/2 
and m^^^ and m™^^) are bounded from below by m|'"'^ (m^™°), while m™^^ (mp*^^) 

is in addition bounded from above by £'miss/2 To compare the difference between 

these two methods over a range of masses, we show two additional SUSY models, specified by 
(mf™,m|"^) = (100,50) GeV and (m*™,m|"^) = (200, 100) GeV in addition to the 'reference' 
SUSY model and the SM background. We see that, compared with the {E^iss^m^^) distribu- 
tions, the (m™^^,m?"^^) distributions benefit from a better signal-to-background ratio (S/B) in 
the signal region, which improves as the mass scale increases. Here, the signal region for a distri- 
bution is defined as the region above the lower threshold value determined by (m^''™, m|^"*^). For 
example, for (m^''™, m|'"*^) = (150, 100) GeV, the S/B ratio in the signal region is approximately 
1 for both TTi™^^ and m™*^^ distributions, while it is approximately 1/4 and 1/3 for the -Emiss and 
m^^ distributions respectively. Importantly, the distribution of the signal -Emiss distributions 
near the endpoint becomes increasingly smooth as the mass scale increases, leading to possibly 
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(left) and m^^ (right) distributions. No detector smear- 



large uncertainties in determining the endpoint value. In contrast, the impact on the m^^^ 
distribution is relatively mild. 

To estimate the precision in determining (mi''"'^, m|'"'^) that can be obtained by measuring 
^^max^ j^max^ LHC, we perform a pseudo-experiment which includes detector smearing 

effects. The 4-momenta of the final state leptons are smeared via 

^ = ^ + ''% (20) 
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where a = 0.1 vGeV and b = 0.007 45|, |46|. The fractional momentum loss C, of the final state 
protons is smeared via0 

J = 0.04. (21) 

In the pseudo-experiment, the number of events we generate for both the SUSY signal, 
^j^true^ j^toue-j _ (^150^100) GeV, and SM background is based on the cross sections displayed 

in Table [H assuming an integrated luminosity of 300 fb~^, i.e. 216 and 38 events for the SM 
and SUSY model respectively. We further define the signal region as the range within m™^^ = 
[130 : 230] GeV and m™*^ = [80 : 180] GeV to increase the S/B ratio. In our particular pseudo- 
experiment sample, 24 and 36 events are observed in the signal region for the SM and SUSY 
model respectively. 

We construct the posterior probability density of (mr, m^) and identify the 68% and 95% 



credibility region as follows. Given a hypothesis m = (m^-, m^^), the probability density function 
(p.d.f.) for the SUSY sample in obtaining m™^'' = (m™^"^, m™'^''), denoted by psusY(m™'^''; m), 
is estimated by generating 1 million events, and binning the events passing the cuts into a grid 
with bin size 2 GeV x 2 GeV on the (m™'^^, m™'^^) plane. We generate 1189 SUSY mass points 
in the range of m^- = [145 : 152] GeV and = [93 : 103] GeV with a 0.25 GeV step size. The 
p.d.f. for the SM background sample, psM(m™'^^), is estimated by generating 1 million events in 
the same way. The combined probability density, p(m'""^; m), is then given by 

p(m--^;m) = /sMPsM(m--^) + /susYPsusY(m--^; m) , (22) 

where fsu and /busy are the fraction of SM and SUSY cross section after the cuts, with /sm + 
/busy = 1- 

The likelihood function, i.e. joint p.d.f., may be written as 

L(m) = JJp(m"^^";m) , (23) 

i 

where i runs over the number of events observed in the signal region. We define the posterior 
probability density a^ 

, , L(m) , , 

P(m) = ... , (24) 
J dm' L(m') 

A hypothesis m is inside the 68% (95%) credibility region if 

j dm!p{m!) < 0.68 (0.95), (25) 

where the integration is carried out over m' satisfying p{m') > p{m). 

^We choose this level of smearing to reproduce the quoted 'missing mass' resolution of roughly 2 — 3 GeV for 
a 120 GeV SM Higgs boson Q 

^ In other words, we take a constant prior p.d.f., 7r(m) = c, in the Bayesian posterior p.d.f. 

L(m)7r(m) 
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/ dm'L(m')7r(m') 
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Figure 6: Credibility distribution for the {mj, rrij^) hypothesis in a pseudo-experiment. The mass 
parameters of the SUSY model are given by (m^™*^, m|'"*^) = (150,100) GeV. The solid black 
and dashed red contours contain the 68% and 95% credibility regions respectively. 



We show the p(m)Am;-Am^ distribution (Am^ = Am^ = 0.25 GeV) from an event sample 
from one pseudo-experiment in Fig. [6l In this particular sample, with 60 events (24 SM back- 
ground, 36 SUSY signal) passing the selection criteria, (m^''"°, m|'"'^) can be estimated to about 
3-4 GeV accuracy at 95% credibility. Results from other event samples, which assume the same 
integrated luminosity, show similar level of accuracies. By applying the likelihood calculation on 
event samples which do not include finite detector resolution effects, we find that the uncertainty 
is primarily due to statistical fluctuations. 

As can be seen in Fig. [6l there is a preference towards the lower mass hypothesis. This is a 
consequence of the peaking structure of the {mj^^^, m™^^) distributions just above {m}- 



This implies that a hypothesis (m[, m^) where m^- > mj' 



true 



many signal events in the region where psusy(™ 
likelihood contribution from /smPsm("i? 



and rriy. > m|'"^ 



true ^true^i 



is likely to have 



max ,™ max > 

I ' '"x 



likelihood reduces less quickly in the case where m^- < rrij^^^ 



; m) is zero, and therefore they only 'see' the 
which is rather low. On the other hand the 
and m 



the mf^^^ and m^^"^ 



distributions. 



< m^''™ due to the tail of 

A. 



4 Summary and discussion 

In this paper, we have discussed how the additional kinematic information obtained in CEP 
with tagged forward protons can be used to simultaneously measure the masses of the particles 
involved in the short decay chain displayed in Eq. ([5]). We presented a full solution to the system, 
and found that the distributions of the quantities m™^^ and m™*^^, which are the maximum mass 
values consistent with a given event, are particularly useful in estimating their corresponding 
true mass values m~^^^ and ml'™. We then performed a likelihood analysis to show that for our 
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reference SUSY model, (m^''"°, m|'"'') = (150, 100) GeV, the true masses can be estimated to a 
few GeV at 95% credibility with just (9(30) signal events. 

As emphasised in Ref. |35|, the possibility to measure central exclusive production in the high 
luminosity phase of the LHC is crucial to a wide range of processes studied in the literature. 
However, one potentially important issue which our simulation has not included is the effect 
of pile-up: at high luminosity there are expected to be multiple overlap events in each bunch 
crossing which will fill the rapidity gaps between the central system and the outgoing protons. 
Even so, by making use of the forward proton taggers and the high precision time-of-flight 
detectors discussed in Ref. |28|] to demand that the tagged protons and the opposite-sign lepton 
pair all point to the same origin, we may hope to suppress the effect of pile-up, given good 
vertexing resolution. However a more detailed investigation is beyond the scope of the present 
study. 

Nevertheless, we find the level of accuracies obtainable with the small (0(30)) event sample 
encouraging. The prospects for measuring the masses in the corresponding decay chain in KK 
scenarios are even better, due to the much higher signal cross section given the same mass 
parameters, see Fig. [TJ 

As discussed in the Introduction, the production of a slepton pair, followed by the decay 
into neutralinos and leptons in an inelastic scattering is clearly possible. As is well-known, the 
energy and the longitudinal momentum of the hard process cannot be measured in a hadron- 
hadron collision. Compared to the CEP process that we consider, the inelastic process has a 
much larger cross section. In principle, the MT2 method can be used to determine the true 
masses by observing a kink structure at the true mass value. This kink structure is formed by 
events where the two slepton system recoils against initial state radiation and has a non-zero 
transverse momentum. However the probability that the system has a significant transverse 
momentum is very low and the kink resides at the smooth tail of the distribution. Once detector 
smearing effects are included, it is likely that observing such a kink structure will be difficult. In 
this context, measuring the particle masses in this particular decay chain in CEP with forward 
proton tagging might be the best, if not the only, way. 

Finally, we note that the mass measurement technique presented here could also be used in 
future linear colliders. A relevant process could be 



e 



^ Itl 



1 '■2 1 



(26) 



where the p^^ is simply the vector sum of the incoming electron 4-momenta. At parton level, 
the leptons have a flat energy distribution with endpoints -Emax.mm at 



max, mm 



4 



true 



1 ± 



2^truc 
7^ 



(27) 



cf. Ref. [50[. Since the masses are determined by the endpoints, the events away from the 



endpoints are not directly used. Our m^^^ and m™^^ distributions on the other hand peak 
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strongly near m^-^^'^ and m|^"°, so they are likely to provide better statistics given the same 
number of events. However we postpone a comparative study of these two measurement methods 
to the future. 
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A Differential QED cross sections with 77 initial states 



For ease of reference, the differential cross sections for pair production of final state particles for 
1^+1^-0, ^+^-13 and are listed below: 



da 



dcos 9* 
da 



dcos 6' 
da 



(77 ^ W+W-) 

ill -> s^s-) 



TTa^f3 19 - 6/32(1 - /32) + 2(8 - 3f3^)(5'^cos'^ 6* + 3f3\os'^ 9* 



(1 -/32cos2 9*y 
2(1-/32) 2(1-/32)2 



(1 -/32C0S2^*) (1 -/32C0S2^*)2 



dcos 9* 



(77 F+F' 



2na^l3 1 + 2/32(1 - /32)(1 - cos^ 9*) - P^cos^ 9* 
~s (1 -/32cos2^*)2 



f28) 



where /3 = {l—Am^/sY^'^, (m = mw, rnp, nis), 9* is the angle of out-going particles with respect 
to the 7's in the CM. frame. The charges and colour factors are all assumed to be 1. The total 
cross sections are given by 



cr(77 ^ W+W- 
a(77 ^ S+S- 
(t(77 ^ F+F' 



7ra2/3 



27ra2/3 



47ra2/3 



■^22-9/32 + 3/34 A-I3\ fl + f3 
2 — 3 — - — In 



1-/32 



2-/32 



1-/34 
2/3 



In 



/3 
1 + /3 
1-/3 



1-/3 



-2 + /32 + ^lnri±^ 

' 2/3 \l-(3 



(29) 
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